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Momentum conservation in dissipationless reduced-fluid dynamics 
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The momentum conservation law for general dissipationless reduced-fluid (e.g., gyrofluid) models 
is derived by Noether method from a variational principle. The reduced-fluid momentum density and 
the reduced-fluid canonical momentum-stress tensor both exhibit polarization and magnetization 
effects as well as an internal torque associated with dynamical reduction. As an application, we 
derive an explicit gyrofluid toroidal angular-momentum conservation law for axisymmetric toroidal 
magnetized plasmas. 



Nonlinear reduced-fluid models play an important role 
in our understanding of the complex dynamical behav- 
ior of strongly magnetized plasmas. These nonlinear 
reduced-fluid models, in which fast time scales (such as 
the compressional Alfven time scale) have been asymp- 
totically removed, include the reduced magnetohydrody- 
namic equations [lH3|, the nonlinear gyrofluid equations 
d, Q , and several truncated reduced-fluid models (such 
as the Hasegawa-Mima equation 6, 7] and the Hasegawa- 
Wakatani equations [8]). Because the space-time-scale 
orderings for these reduced-fluid models are compatible 
with the nonlinear gyrokinetic space-time-scale orderings 
Q , they provide a very useful complementary set of equa- 
tions that yield simpler interpretations of low-frequency 
turbulent plasma dynamics in realistic geometries. 

The self-regulation of anomalous transport processes 
by plasma flows in turbulent axisymmetric magnetized 
plasmas has been intensively investigated in the past 
decade. Because a strong coupling has been observed [|[ 
between toroidal-momentum transport and energy trans- 
port in such plasmas, it is natural to investigate the 
link between these two global conservation laws through 
an application of the Noether method on a suitable La- 
grangian density [Io| . The purpose of the present Letter 
is to focus its attention on a momentum conservation 
law derived from a general reduced-fluid model [llj and 
then explicitly investigate the reduced toroidal angular- 
momentum transport in axisymmetric magnetic geome- 
try derived from it. 

The general variational formulation of nonlinear dissi- 
pationless reduced-fluid models is expressed in terms of 
a Lagrangian density C(ip a ) as a function of the multi- 
component field 



i> a = ($,A,E,B;n,u,pi|,px) 



(1) 



Here, the electromagnetic fields (E, B) are defined in 
terms of the electromagnetic potentials ($, A) as 

E = — V$ — c^dA/dt and B = V X A (2) 

and the reduced-fluid moments (n, u,pu,px) are used for 
each plasma-particle species (with mass m and charge 
q). We note that the Lagrangian formalism does not ac- 
commodate higher-order fluid moments (e.g., heat fluxes) 



and, therefore, the issue of fluid closure is completely 
ignored [12]. These higher-order moments, as well as 
dissipative effects, can be added after the dissipation- 
less reduced-fluid equations are derived by variational 
method [l3[ (although a variational procedure [HI, HH 
can be used to include heat fluxes in the pressure evolu- 
tion equations). 

We begin with the general reduced Lagrangian density 



C(iJ; a ) = £ M (E,B) + £*($,A;n,u) 
+ £p(n,u,p||,px;E,B), 



(3) 



where the electromagnetic Lagrangian density is Cm = 
(|E| 2 — |B| 2 )/87r, the gauge-dependent interaction La- 
grangian density (summed over particle species) is = 
— 1 n — A • u/c), and the reduced-fluid Lagrangian 
density £f depends on (E.B) only through the process 
of dynamical reduction [16| . 

The reduced plasma-electrodynamic equations associ- 
ated with the reduced Lagrangian density ([3]) are divided 
into either constraint equations or dynamical equations. 
The electromagnetic constraint equations are 



V-B 







V X E 



^dB/dt, 



(4) 



which are satisfied by the representation @. The 
reduced-fluid constraint equations, on the other hand, 
are the continuity equation 



dn „ 



n u 



(5) 



and the Chew-Goldberger-Low (CGL) pressure equations 

mm 



= - V- pi: u 



dp 
~dt 

— = -V.^pxu 

associated with the CGL pressure tensor 



- 2p\\ b b : Vu, 

- p ± (I - b b ) : Vu, 



P\\ b b 



P± 



bo b ) 



(6) 
(7) 

(8) 



where Bo = Bq bo denotes the quasi-static background 
magnetic field. The reduced-fluid velocity u appearing 



2 



in Eqs. ©-([7]) will be determined from the variational 
principle 



SCd 4 x = 0. 



(9) 



The reduced Maxwell equations and the reduced-fluid 
momentum equation are the dynamical equations derived 
from the reduced variational principle ([9]). First, the re- 
duced Maxwell equations 



V X H 



V • D = 4tt q, 
1 <9D 
c ~dt 



47T 

"J ■ 

c 



(10) 
(11) 



are expressed in terms of the reduced charge density 
q = — dCq, 1 9$ = 1 n an( i the reduced current den- 
sity J = cdC^/dA = qnu, while the reduced elec- 
tromagnetic fields D = 47r<9£/<9E = E + 4irP and 
H = — Att dC I <9B = B — 4ir M are expressed in terms 
of the reduced polarization and magnetization 



(P, M) 



{ <9E ' dB 



(12) 



Equations (fT0|) -([TT |) can also be expressed as 

V-E = 4tt (q - V-P) , (13) 

VxB-i£-£( J + £ + cVxM), ( U, 

where g po \ = — V • P denotes the polarization density, 
Jpoi = dP/dt denotes the polarization current, and 
Jmag = cV X M denotes the magnetization current. 
Next, the reduced-fluid momentum equation is 



dp 

~dt 



qn E 



xB 



n{VK - Vu-p) 
+ V-P.V (15) 



- [P± V 7 ± + — V 7 || 

where d/dt = d / dt + u • V , the reduced-fluid kinetic en- 
ergy K and the reduced-fluid kinetic momentum p are 



K 
P 



dC F /dn 
n^ 1 dC-p/dn 



(16) 



and the symmetric reduced pressure tensor 

P* = p||7||b bo + p ±7± (I-b b ) (17) 

is defined in terms of the coefficients 7|| = — 2dC-p/dp\\ 
and 7^ = — dC-p / dp± . We note that this pressure ten- 
sor generalizes the CGL pressure tensor @ and includes 
standard finite-Larmor-radius (FLR) corrections through 
7X^1 El. 

The reduced equations Q-Q, ([TO ]) - ([IT ]) . an dfPo]) sat- 
isfy the reduced momentum conservation law [11| 



on 

~dt 



V-T 



V'£, 



where the reduced momentum density is 

D x B 



n 



47T C 



(19) 



the reduced canonical momentum-stress tensor is 



T = H P * 



IEI 



1 + \B\ 2 ) - B-M 



8tt 

— (e E + B B 

47T V 

[^nup - fPE BM 



(20) 



and V'£ denotes the spatial gradient of the reduced La- 
grangian density £ = £ — with the dynamical fields 
((T|) held constant. We note that, while the first three 
terms in the canonical momentum-stress tensor (|2"0)l arc 
symmetric while the remaining terms (on the last line) 
are not. The antisymmetric part (Ta)^ = \ (Ty — Tji) = 
of the canonical momentum-stress tensor (1201) 



can be expressed in terms of the reduced internal torque 
density 



t = ^ nuxp + (ExP IBxM 



(21) 



which exhibits classical zitterbewegung effects jl Qj j associ- 
ated with the decoupling of the reduced-fluid momentum 
p ^ mu and the reduced-fluid velocity u and the reduced 
polarization and magnetization effects. Here, the inter- 
nal degrees of freedom of a gyrofluid particle are associ- 
ated with the fast gyromotion that has been eliminated 
by dynamical reduction. 

The symmetry of the momentum-stress tensor is phys- 
ically connected to the conservation of angular momen- 
tum [l(| H3, i.e., conservation of the total angular mo- 
mentum (including internal angular momentum) explic- 
itly requires a symmetric momentum-stress tensor. Since 
the left side of Eq. (TT81 is invariant under the transfor- 
mation [H| n' = n + V-S and T' = T - dS/dt, the 
second-rank antisymmetric tensor S = | Sijk <J k can be 
chosen so that T' = T$ = \(T%j + Tji) is symmetric, i.e., 
the antisymmetric part Ta = dS/dt yields the reduced 
internal angular momentum equation 



da 

~dt 



(22) 



where <x denotes the internal (spin) angular momentum 
density, and the reduced momentum conservation law 
(11811 becomes 



^(n-^VXff) f V-T s = V'£. (2:-!) 



(1^) where we used the identity V • S = — iVxer. 
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By applying the Noether Theorem in axisymmetric 
tokamak geometry, where 



Bo = V<pxVtp + B 0ip (iP)V<p 



(24) 



and the background scalar fields are independent of the 
toroidal angle ip (i.e., d'C/dip = 0), we obtain the reduced 
toroidal angular-momentum transport equation fill ] 



d i 



Here, the toroidal momentum density is 



0. 



(25) 



n ^ x 

dip 



z-x x n 

47T C 



(26) 



D X B 

47T C 



<9x 

dp' 

fell Bn 



where the perturbed magnetic field is B — Bo 
B^ and we used the identity B X dx/dp = Vip. 

A more useful expression for Eq. (|25[) . however, is ob- 
tained in terms of the magnetic-flux average (• • • ) = 
V- 1 § J (■■■) d9 d<p as 



dt 



- — (V (t* ) 



(tz), (27) 



where J = (Vtp X V6> • Vy>) _1 = (B • VA)" 1 is the Ja- 
cobian for the magnetic coordinates (ip,0,ip) and V = 
§ J dd dp. In Eq. (|2T|) . (r z ) denotes the reduced in- 
ternal torque density and the surface-averaged toroidal 
angular-momentum flux is 

9x 



r 



ip 



(28) 



1 

47T 



where B]_ = Bi • denotes the V'-component of the 
perpendicular component of the perturbed magnetic field 
(since B • W> = 0). 

We now investigate the surface-averaged toroidal 
angular- momentum conservation law (1271) by considering 
the gyrofluid Lagrangian density [H, |22| 



E 
E 



IBI 



— mn tin 
2 11 



(nK p + V) 



qn 



c V 



■ bo 



(29) 



where uii = u • bo denotes the gyrofluid velocity along the 
unperturbed (background) magnetic-field lines and the 
pressure tensor (fTT|) is P* = P (i.e., 711 = 1 = 7_l) with 
V = \ Tr(P) = p x +p\\/2. In Eq. {25]), the zero-Larmor- 
radius limit of the gyrocenter dynamical reduction [f| 
introduces the nonlinear FLR-corrected potentials 



A,, 




(30) 



which generate the nonlinear FLR-corrected electromag- 
netic fields 



B, 



- W p - c- x b dA llp /dt 
V x b ) 



and the low-frequency ponderomotive potential 



(31) 



(32) 



which generates the generalized ponderomotive force 
density V • P p = V • P + n V/C p . These nonlinear FLR- 
corrected fields are expressed in terms of the gyrofluid 
displacement [22| 



P± 



Ei 



b X B j 



bo 



X U, 



(33) 

The partial derivatives (| 16 j) of the gyrofluid Lagrangian 
density (|29|) yield the gyrofluid kinetic energy K — 
m + PCp, and the gyrofluid momentum 

p = m + U ± -^-j b = muj b , (34) 

where itjj defines the gyrofluid velocity along the per- 
turbed magnetic-field lines. Lastly, the partial deriva- 
tives (|12l) yield the gyrofluid polarization and magneti- 
zation 



P = ^2 qn p ± = E TOn 7T x U± ' 



Bo 



(35) 



qn p , X — b = > mn -f Ui, (36) 
c ^— ' Bq 

which appear in the Maxwell equations (|10p -(fTT j) . Note 
that the reduced internal torque (1211) for this gyrofluid 
model is expressed as 



= E 



(uy — ut) I E 



'II 



bo X B± 



mn ut ( Ui — j X b 



(37) 



where E_l X P + B ± X M = 0. 

The gyrofluid momentum equation (|15l) derived from 
the gyrofluid Lagrangian density (|29p is expressed as 

mn b = qn (e p + " xB p ) - V-P p , (38) 

where B* = B + U|| (B /Q ) V X b + B ip . Equation 
(|38| contains both the gyrofluid parallel-force equation 



dun 
~dt 



b* p ■ (qn E p - V • P, 



where b* = B*/£> , and the gyrofluid velocity 



mii b* 



(qn Ej_ p 



mn i Iq 



(39) 



(40) 
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which includes the E x B velocity and the diamagnetic 
velocity and their nonlinear FLR corrections. 

Next, we consider the gyrofluid version of the surface- 
averaged toroidal angular-momentum conservation law 
([27]) . where the gyrofluid momentum (l34l) is substituted 
in Eqs. (p ]) -([20 ]l . with P» = P and r] a dC F /dTj a = 
£f in Eq. (j2TJ)) . The gyrofluid version of the surface- 
averaged equation (|27p is expressed in terms of the gy- 
rofluid toroidal angular-momentum density (|26)). with 



mU|J bo v and 



Emnc 2 x 



R2 
-°0 



mnu||C N 



R2 



B 



.-(V^Xbo), (41) 



while the expression for Eq. 



n - e 



qp]_E v + mun — U± v 



4?r 



B^B V J 



(42) 



where we have separated the reduced polarization P^ = 
1 n P± and magnetization — mnu\\U ^ / Bq 

from the Maxwell stress tensor. 

As an application of the gyrofluid model (f2U|) . we con- 
sider its electrostatic version (E = — V<f>, 



II 



'II 



and B = Bq) and use the quasi-neutrality condition 



g = V • P [valid for An (J2 mnc 2 / B%) > 1]. The surface- 
averaged gyrofluid toroidal angular-momentum equation 
(1271) therefore becomes 



5 (<Hh> + !<"♦>- W) 



v 



pi/; 



dip 



(43) 



where = mnu||)& 0y , = (J2 mnu^u t( ')b 0v , 
and = (J2 mnc 2 / Bq) E^ . This equation was re- 
cently obtained [23[ (without the internal angular mo- 
mentum <j z ) by direct evaluation of the time evolution 
of the surface-averaged gyrocenter moment (n^ an ) = 
(n^ii) — (jjj/c) (g) of the toroidal canonical momentum 
mv \\ boip — qip/c, where the surface-averaged gyrofluid 
charge density (g) = V -1 d(V (P^)) / dtp is expressed 
in terms of the surface-averaged polarization component 
{F*). 

In this Letter, we have derived an exact toroidal 
angular-momentum conservation law (|25p that clearly 
exhibits the role played by the reduced internal torque 
t z = da z /dt in possibly driving spontaneous toroidal ro- 
tation in axisymmetric tokamak plasmas. 
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